Nonadditive Quantum Error Correcting Codes Adapted to the Amplitude Damping 

Channel 
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A family of high rate quantum error correcting codes adapted to the amplitude damping channel 
is presented. These codes are nonadditive and exploit self-complementary structure to correct all 
first-order errors. Their rates can be higher than 1/2. The recovery operations of these codes can 
be generated by a simple algorithm and have a projection nature, which makes them potentially 
easy to implement. 



1. INTRODUCTION 

As quantum computation finds wide applications to- 
day, the difficulty of implementing quantum computers 
also arises from the expense of qubits and their vulnera- 
bility to dccohcrcncc. Therefore, quantum error correc- 
tion and fault tolerantquantum computation have been 
extensively developed [9j . Among the most famous quan- 
tum error correcting codes arc the [9,11 code H3, the CSS 
[7,1] code 0,03, and the [5,1] code 0,0]. The encoding 
rate of these codes are much lower than classical codes 
due to the complicated nature of quantum decoherence. 
However, qubits are expensive, so it is desirable to have 
codes with higher rate. Unfortunately, it has been shown 
that 1/5 (encoding 1 qubit into 5 qubits) is the highest 
rate for one-qubit code. In fact, quantum error correction 
codes satisfy the following quantum Hamming bound Q : 
if an [to, k] code (encoding k qubits into n qubits) can 
correct errors at at most t qubits and there are a inde- 
pendent possible errors in one qubit, then 
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(1.1) 



It has been shown that if a code corrects all the errors in 
the Pauli group generated by the Pauli matrices operat- 
ing on each qubit, then it can correct all types of errors 
0. Therefore, traditional quantum error correction has 
concentrated on the errors in the Pauli group, and hence 
a = 3 in (jl.ip . Within this framework, it is hard to 
improve the rate further. 

However, not all the errors are equally likely in a re- 
alistic channel. If we assume some knowledge about the 
channel, such that only 1 or 2 independent errors are 
likely to occur in each qubit, then codes of higher rate can 
be obtained. Indeed, an iterative optimization method to 
investigate channel-adapted codes has been proposed by 
[J] and Q and a family of [2n + 2, to] stabilizer codes, 
which correct all the first-order errors and some of the 



second-order ones in an amplitude damping channel, have 
been discovered by The amplitude damping channel 
is the tensor product of one-qubit amplitude damping 
channels, each of which consists of the following two op- 
erator elements: 
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(1.2) 



where 7 is a parameter measuring the strength of the 
noise. Here |1) is pictured as an excited state and |0) the 
ground state. 7 is the probability of the transition from 
|l)to|0). _ 

Another important concept .nonadditive codes, was in- 
troduced in recent research 0, [HI, [13] seeking to find 
codes with higher rates than stabilizer codes. In these 
codes, the codewords do not form a subspace of F£ and 
thus the dimension of the source space is not of the form 
2 . In other words, these codes encode a fractional num- 
ber of qubits. To avoid confusion, we denote a (nonad- 
ditive) code which encodes a k dimensional space into n 
qubits by an (to, k) code. This paper combines the ideas 
of channel adaption and nonadditive codes to find a new 
family of codes adapted to the amplitude damping chan- 
nel which outperform the [2to + 2, to] stabilizer codes in 



2. SELF-COMPLEMENTARY CODES 

The family of codes to be presented exploits self- 
complementary structure, which is crucial for correct- 
ing errors in the amplitude damping channel. Self- 
complementary codes appear in different contexts 0, [ll| 
and play different roles. To understand the significance 
of the structure, we analyze how errors arise and are cor- 
rected as follows. 

Let C be a quantum channel. The extent to which the 
channel preserves information is measured by the entan- 
glement fidelity [§] 



F(p,C(p))= try/ ptC(p) f 



(2.1) 
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where p is the input state. Sometimes, we are more inter- 
ested in the worst case performance, which corresponds 
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to the minimum entanglement fidelity 

F min (C) = minF(|V) (VI ,C(|$ «>|)). (2.2) 

The restriction to pure states is justified in [9(. In our 
case, C is composed of the encoding operation U, the am- 
plitude damping channel £(7) and the recovery operation 
7^(7). The dependence of 1Z on 7 suggests that the code 
be channel-adapted. The errors contained in C(j), ap- 
propriately called the residual errors, can be divided into 
two classes. First of all, an error may flip a codeword 
\wi) (wi\ to another codeword \w2) (w2\, like an X op- 
erator acting on this two-dimensional subspace. We call 
this an X-like error. Assuming no X-like errors (up to 
a specific order), we may still have an error that causes 
damping or a phase shift (or both) of an off-diagonal en- 
try: 1 — tr(|w2) (wi\C(\wi) (W2D) 7^ 0. Different errors 
cause this problem, including the Z operator in the two 
dimensional space spanned by \w\) and \w2). For sim- 
plicity, we call these errors Z-like errors. The following 
proposition shows that this list is exhaustive, up to the 
first order. 

Proposition 2.3. The code corrects all the first-order 
errors in the sense that F m - m = 1 — o(7 2 ) if and only if 
£(7) do not contain any X-like errors or Z-like errors to 
the first order. 

Proof. If £(7) contains an X-like error on \w±) or a 
Z-like error on (102 |j then choose — \wi) or 
\ip) = {\uii) + |u; 2 ))/v / 2, respectively, to show that 
F(\ip) (ip\ ,C(\tp) (ip\)) contains a first-order term. 

Conversely, assume that £(7) does not contain any X- 
like or Z-like errors to the first order. Call the underlying 
Hilbert space H, which is finite dimensional. Then C(j) 
is a linear operator on C(H,H), the space of the linear 
operators on H. By the assumption, the entries of the 
matrix / — C(j) in the basis {\wi) (1U2I} do not contain 
first-order terms, so j|/ — C (t) 1 1 = o(7 2 )- The conclusion 
follows from the continuity of F(p, a) and the compact- 
ness of the unit sphere in H . □ 

Consider an (n, k) code. For a codeword u G {0, 1}™, 
we denote by u its complement such that u j = 1 — Uj for 
all j. We call a code self-complementary, if the code is 
spanned by kets of the form (|u) + \u))/V2. We define 
the inner product of two codewords u and v by u ■ v = 
Sj=i u j v j an( i thus the Hamming weight of u is \\u\\ 2 . 
By the definition (|1.2j) of the amplitude damping channel, 
after passing through the channel, (\u) + \u))/\/2 becomes 
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|u) + (l- 7 )i 



|«)). (2.4) 



Let \f{u)) be the unit vector associated to f(u)), and 

\g(u)) be the unit vector in span{|u) , |u)} orthogonal to 
|/(it)). It is thus reasonable to correct \f(u)) to the 
codeword (|ti) + \U})/\/2. When we talk about self- 
complementary codes adapted to the amplitude-damping 



channels, wc will always assume this recovery operation. 
Now we will show that the self-complementary structure 
eliminates the Z-like errors provided that the X-like er- 
rors have been eliminated. 

Theorem 2.5. Let C be the composite channel of a 
self- complementary code subject to an amplitude damp- 
ing channel. If £(7) does not contain X-like errors to 
the first order, it does not contain Z-like errors either. 

Proof. Let e«' € {0, 1}™ denote the error vector such 



that ef = Sij. By the assumption, if \w 
\u))/y/2 is in the code, then |u 
fully to \w) provided that Uj ■ 
results in a first-order X-like error. When 



ev)\ will be corrected 
1 since otherwise this 



agree that 



) 



= 0, we 

0. Now consider two codewords 



K) = (|u) + \u))/V2 and \w 2 ) 
u, v G {0, 1}™. We have 



(\v) + \v))/y/2, where 



C(\w 1 )(w 2 \) 



U) 



K[f(u)){f(v) +lJ2(\u-e 
+ u - e {: >A (v - e (j) +\u- e ij) 



) + o{l 2 )] 



v — e 



u — e^' ) < v — 



^((1 _ 7 )||«p + (i _ 7 )INI a )((l - 7)NP + (1 - ipV) 



x \w\) (ti^l + — (u ■ v + u ■ v + u ■ v + u ■ v) \w\) (102] + o(7 2 ) 

= i(2-^(|| U || 2 + ||^|| 2 + || W || 2 + ||5|| 2 ))| U ; 1 )( W 2| 
+ ^(w + u) ■ (v + v) \wi) (w 2 \ + o(7 2 ) 
Since ||m|| 2 + ||u|| 2 = ||u|| 2 + ||i;|| 2 = (u + u) ■ (v + v) = n, 



the conclusion follows. 



□ 



CODEWORDS AND RECOVERY 



It is clear from Theorem 12 .51 that a self-complementary 
code corrects all the first-order errors in an amplitude- 
damping channel if and only if no confusion arises assum- 
ing the decay occurs at no more than one qubit. More 
precisely, C = span{(|u) + \u))/\/2\u € S} is such a code 
if and only if the set S C {0, 1}" satisfies 

(51) If u e S, then u e S; and 

(52) If u, v G S and u — = v — e^> for some then 

u = v. 

The dimension of the code is k = dimC = |<S'|/2. We use 
a greedy algorithm to search for maximal sets S satis- 
fied the above conditions for some small values of n and 
the results are listed in Table 1 . The log 2 k column ap- 
proximately indicates the number of qubits that can be 
encoded. Linear regression yields that the slope of y ~ x 
curve is 0.85, which is higher than 0.5 as in [5j. 
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Tabic 1: Encoded Dimensions for Small Values of n 



n 


k 


1 7 

log 2 k 


4 


2 


1.00 


5 


2 


1.00 


6 


5 


2.32 


7 


8 


3.00 


8 


12 


3.58 


9 


18 


4.17 


10 


41 


5.36 


11 


78 


D.ZO 


12 


146 


7.19 


13 


273 


8.09 


14 


515 


9.01 


15 


931 


9.86 


16 1716 


10.74 



We now describe the algorithm to generate the recov- 
ery operation. We use the maximum likelihood recovery. 

Algorithm 3.1. The error vectors refer to vectors in 
{0, 1}™, where 1 indicates a decay in that qubit. The 
error vectors are sorted by their weights and those with 
the same weight are sorted in the dictionary order. 

Step 1 The first operator element of TZ corrects | f(u)) 
in (|2.4[) to (| it) + \u))/^/2. We are considering the 
error vector e = (0, 0). 

Step 2 Put e to be its successor and come to construct 
the next operator element oflZ. 

Step 3 Find a word u € S such as u — e > 0, which 
means that the decay e may happen to u, do the 
following: 

Step 3.1 If u — e e S and both u — e and u — e 
have not appeared in TZ except in the first 
operator element, then correct \g(u — e)) to 
(\u) + \u))/y/2. Go to Step 3.3. 

Step 3.2 If u—e has not appeared in TZ, then correct 
\u-e) to O) + \u))/y/2. 

Step 3. 3 Search for the next word u € S with u — 
e > and go back to Step 3.1, until all the 
words in S have been exhausted. 

Step 4 If the sum of the ranks of all the constructed oper- 
ator elements ofTZ equals n, then stop. Otherwise, 
go back to Step 2. 

By the construction, the operator elements of TZ are 
orthogonal projections, which makes it easy to implement 
as a quantum circuit. By Theorem l2.5l TZ corrects all the 
first-order errors in the amplitude damping channel. 



The Entanglement Fidelity of the (8, 12) Code 
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FIG. 1: The entanglement fidelity of the (8, 12) code versus 
the parameter 7, on the initial state po = 1/12. The dashed 
curve is the entanglement fidelity of the bare 3 qubits on the 
initial state 1/8. 



4. AN EXAMPLE 

We demonstrate the performance of our codes by a 
simple example, the (8,12) code, which encodes a 12 di- 
mensional space to 8 qubits. In this case 

5 = {00000000,00000011,00001100,00110000.11000000, 
10101000, 01011000, 01100100, 10010100, 11110000, 
11001100, 00111100 and their complements} 

The entanglement fidelity of this code on the ensem- 
ble J/12 is plotted and compared with the unprotected 
3 qubits, since this code can encode at most 3 qubits. 
We can see that the code has the desired behavior when 
7^0 (the linear term vanishes) and it also has a good 
performance for larger 7. See Fig. [T] 

5. CONCLUSION 

We have described the structure of a family of non- 
additive quantum error correcting codes adapted to the 
amplitude damping channel with the emphasis on the sig- 
nificance of its self-complementary structure. The code 
has an extremely high rate and thus may turn out use- 
ful in quantum computation and quantum teleportation. 
Since the code is nonadditive, its decoding docs not have 
the syndrome-diagnosis and recovery structure of the sta- 
bilizer codes, so we are concerned with designing the effi- 
cient quantum decoding circuits for these codes. It is also 
challenging to figure out the cardinality of S satisfying 
(SI) and (S2) mentioned in the third section in a closed 
form or its asymptotic behavior. Table 1 is constructed 
using the greedy algorithm, so S may well have a larger 
size than listed there. 
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